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Abstract 

By allowing interactions between electrons, a new Harper's equation is derived to ex- 
amine the effects of electron correlations on the Hofstadter energy spectra. It is shown 
that the structure of the Hofstadter butterfly for the system of correlated electrons is 
modified only in the band gaps and the band widths, but not in the characteristics of 
self-similarity and the Cantor set. 
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1 Introduction 

Numerous studies have been made to study two dimensional quantum systems of electrons 
in a magnetic field, including the Hofstadter energy spectrum [1] and the quantum Hall 
effects [2]. Hofstadter studied the energy spectrum for non-interacting electrons on a two 
dimensional square lattice in a perpendicular uniform magnetic field, namely the energy 
spectrum vs. magnetic flux per plaquette known as the Hofstadter butterfly. The energy 
spectrum was found to critically depend on the ratio p/q {p and q are positive integers) of 
the magnetic flux per plaquette to the flux quantum. If p/g is a rational number, each energy 
band is split into q sub-bands by the magnetic field. The Hofstadter butterfly displays a 
recursive structure over rational (number) fields and a Cantor set at any irrational number 
of p/g [1]. 

The gap structure predicted by the Hofstadter butterfly is known to play an important 
role on the magnetoresistance of the laterally modulated two-dimensional electron systems 
in GaAs-AlGaAs heterostructure [3]. Such gap structure can be sensitive to the strength of 
interactions between electrons. Lately, based on the Hartree approximation Gudmundsson 
and Gerhardts [4] studied the effects of interacting electrons on the Hofstadter butterfly. 
In the present study we derive a new Harper's equation which will allow us to study the 
Hofstadter energy spectrum as a function of Coulomb repulsion by properly taking into 
account electron interactions. As an illustration we pay a special attention to the interesting 
system of antiferromagnetically correlated electrons on a square lattice. 
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2 Hofstadter Spectrum of Interacting Electrons 

Earlier Hasegawa et al.[5] studied the Harper's equation [6] and the energy spectrum only 
for spinless non-interacting electrons. In the present study we derive a new Harper's equation 
in order to study the Hofstadter butterfly, that is, the energy spectrum vs. magnetic flux as 
a function of interaction strength (Coulomb repulsion U) between electrons. 

In order to treat interacting electrons in a two dimensional square lattice in a perpendic- 
ular uniform magnetic field, we consider the one-band Hubbard model Hamiltonian, 

^ = - E t^JClC^a + UY: " E ^10,,, (1) 

{ij)a i ia 

with 

tij = -t exp (^-^j- f. A • dlj . (2) 

Here C/^(Cio-) is the creation (annihilation) operator of an electron with spin a at site 
i — {irci iy)- i^j) denotes only the nearest neighbor hopping with lattice spacing a. 0o 
is the unit flux. The choice of the Landau gauge, A = B{0, jxa,0) for a uniform magnetic 
field B leads to the hoping integrals, Uj = 't{j^ii,jy){j:,,jy) — t for the x-direction of hopping 
and tij = t{3x,3yB.){3^,jy) = *exp (^i27rj^^) = iexp (^fSttj^J) for the y-direction of hopping, 
where = Ba^. ^ = | with p and q integers is the number of flux quanta per plaquette 
which is equivalent to the gain of phase by a particle as a result of hopping round the closed 
path of a plaquette. 

We now introduce the staggered magnetization at site 
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with Q = ■^), r = a{ix^iy) and the hole doping rate 5, 

Considering a uniform doping and a uniform staggered magnetization, i.e., 5{i) = 6 and 
m{i) = m respectively, we obtain from (1) the following mean field Hamiltonian for inter- 
acting electrons in the square lattice. 



{ij}a iu 



,m 2 1 



2. 



(3) 



Here the exchange correlation is properly introduced [7] . The third term represents a constant 
energy shift. Ignoring the third term, we rewrite (2) explicitly 



H 



Ua:,Jy-l)cr {Jx,Jy) 

crmU 



Ux,jy)cr 



U 



Ux ,jy)cr 

Omitting the third term which represents the shift of constant energy, the Hamiltonian 
is reduced to 



H - -t J2 ['^h.+ijyW'^U-dyh + ^ij.-i,jyhC' 



(5) 



2 ^Ox,iy)<T'-'O^Jy)<7 



Now we denote \{jx,jy)) as a single (one) particle state in number representation. 



ia,J.)) = |iow.)) = 4...)|o) 
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(6) 



and define 

|(jx±l,jy)) = |l0.±l,i,)) = C'o,±l,j,)C{i.J,)|l(j=.,i«))- 

l(ix,Js,±l)) = |l(i.j,±i)) = '^(j,,j,±i)C'(>,j,)|l(>,i,))- (7) 

We now use the Schrodinger equation, H\il)) = E\iIj) where H is the Hamiltonian (5) and 
the eigenstate in number representation. We then obtain from the use of @ and 



-t 



((j.-l, JyW + ((j.+l, JyM + e-'^^^'^^((j., Jy-im + e^^V^^U,, 3y+im (8) 



777 T J 

-(-1)^^^^^((J., JyW = E{{J., JyW 

Here jx = 1, - ■ ■ ,q acts as a coordinate in the magnetic celL The one- dimensional quasi- 
periodic boundary condition is satisfied; {{jx, jy)\'ip) = ^^^^''"'{{ix + <i,iy)\'ip)- There exists q 
plaquettes per magnetic unit cell. Accordingly the energy spectrum has q sub-bands as a 
result of the applied magnetic field. Here for the range of k vector one magnetic Brillouin 
zone is equivalent to one original Brillouin zone divided by q. 

The above Hamiltonian (5) is invariant under the translation jy jy + 2 for the system 
of antiferromagnetically correlated electrons. Thus using the Bloch theorem, 

{{jx, 3yW = e"'y''^y{{jx, 0)\^) for even Jy, 
and 

{{Jx,jym = e^'«"^«((j., for odd j„ 

we write 

{ij^,Jym = e^'y^^yi^ > {{Jx,0m + ^-^{{Jx,imj (9) 

^^Ka,S^ {{j.. om + {{j., m) , . , j{3x,om-{{3x,im \ 



By defining 



:io.a) 



and 



we rewrite (9) above, 



Using (11), we readily find from (8), 



(10.6) 



(11) 



-t 



u 



mil 
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Tlie expression (12) above leads to the following new Harper's equations which are simply 
the coupled one-dimensional quasiperiodic difference equations. 



~t 



P TflU 

u{jro-l) + u{j^+l) + 2 008(277- j^+kya)u{j^) — 

q j Z 

p . , 1 mil 
v{jx-l) + v{j^+l) + 2cos{2'K-j:a+kya)v{j^) 



v{j,) = Eu{j,) (13) 



The coupled difference equations here arise due to the antiferromagnetically correlated elec- 
trons. 

Defining the reduced energy e — E/t, and the reduced energy gap A = mU/t, we rewrite 

(13) above to explicitly show the following recursion relations, 

p A 
u{j + l) = (-e - 2cos(27r-j + kya))u{j) - u{j-l) - —v{j) 

q z 

p A 
""(i + l) = {^-'^cos{2T:-j + kya))v{j)-v{j-l) + —u{j) 

q Z 



or in matrix form. 

U{] + 1) 

v{j + l) 



V 



t -2coii{2Tx^ i + kyo) -1 



A 

2 



A 

' 2 



v 











e-2cos(27r2j + fcya) -1 
1 



vU-l) 



\ 



(14) 



Since there is no possible confusion below, we omit the subscript x ofjx in the expression ([T^). 



To block-diagonalize the square matrix in ( p!4D above, we introduce the following trans- 
formation, 

^ v'U) ) \-p a J y v{j) J 



a p 



where 



a 



An2 



\ 



Thus the square matrix above is block-diagonalized as 
^ -eo-2cos(27r|j + z/) -1 

















V 








eo - 2cos(27r2j + z/) -1 



where eo = ye^ - 

Using (14) and (15) above, we readily find 



n'(j + l) 
u'U) 



eo — 2 cos(27r2j + /c^a) 



V 



1 ^ I u'ij) ^ 



(15) 



(16.a) 



\ ( 



eo — 2 cos(27r|j + A;ya) 



V 







\v\j-l) ) 



(16.6) 



eo here represents the energy spectrum of non-interacting electrons. The energy dispersion 
relation p| for the system of interacting electrons is readily obtained from (16) above, 



(17) 



with the energy gap, Here we reintroduced the first term which was omitted in (5). As 
seen in ( |T7D above, the energy spectrum e shows the variation of energy gap with both the 
staggered magnetization m and the Coulomb interaction U . However the dispersion relation 
(17) indicates that the self- similarity and homeomorphism to the Cantor set are not affected 
by such electron correlations. 

Now using in (16) the translation matrix 



made of q successive products of 



Aj{e; ky) 



— eQ ~ 2 cos{2tx^ j + kyo) —1 ^ 



V 



we obtain 



1 



Q[e\ kyO) 







' .(1) ] 



(19) 



where g{]) = u'{j) or v'{j). Note that in the expression (|T8|), eo is replaced by e due to 
relation (IT7|). Further we note from the Bloch theorem, 

(20) 







%(!)] 


y 9{q) ) 
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V 9(0) J 



Thus we readily find from (19) and (20), 

det(g(e; fc^,) - e^«*^-"7) = 0. 

One of the eigenvalues of Q is complex conjugate to the other in order to have a real value 
of tr Q. Thus we have 

tr Q{e;ky) — 2cos{qkxa). (21) 

The ky dependency of the trace is additively separable as shown by Buther and Brown [9] , 
and thus we write 

TT 

tr (5(e) = tr Q{e; — ) = 2cos{qkxa) + 2cos{qkya). (22) 
Zq 

Since tr Q{e) is simply the q-th polynomial of e, there exist q sub-bands for a given k. 
Here k is readily obtained for given e. e should satisfy the following condition, 

|trg(e)|<4. (23) 

One can determine from the condition (23) whether a state of a given energy e is possible or 
not. Using the relations (22) and (23) above, we compute the Hofstadter energy spectrum 
as a function of | for the square lattice of antiferromagnetically correlated electrons at half- 
filling, that is, 5 — 0. In principle, the self consistent equations for m and /i can be derived for 
a lattice of infinite size [8]. For the case of finite size lattices we obtain m and //. numerically 
for arbitrary values of p and q. It is important to thoroughly check numerical precision of 
the computed staggered magnetization m in order to accurately evaluate dispersion relation 
(17). This is because m in the dispersion relation is very sensitive to the variation olp/q. 
In order to check numerical accuracy, in Fig. 1 we show the computed results of staggered 

9 



magnetization m as a function of Coulomb repulsion U for the two cases of finite (10 x 10, 
20 X 20 and 30 x 30 square lattices) and infinite size lattices. For both zero magnetic field 
and the non-zero magnetic field (corresponding to half a flux quantum per plaquette), the 
computed results of m with the finite size lattices (denoted as solid or open circles in Fig. 1) 
are nearly the same as the exact values obtained from the use of the infinite size lattice 
(denoted as a solid line in Fig. 1). Encouraged by such numerical accuracy, we computed 
the Hofstadter energy spectra of interacting electrons for various values of U by using the 
finite size lattice of 10 x 10. They are shown in Fig. 2 through 4. We find that the structure 
of the Hofstadter energy spectrum for interacting electrons is greatly modified only in the 
band gaps and the band widths (the bars in the figures stand for the band width) which 
depend on the strength of the Coulomb interaction [/, as well displayed in Figs. 2 through 
4. Indeed this is clearly understood from the expression of the energy dispersion relation 
(17). Although not shown here, at higher values of U we find no other changes but larger 
band gaps and narrower band widths in the Hofstadter energy spectrum. The self-similarity 
of the Hofstadter butterfly is preserved at all values of Coulomb interaction U. 

The band gap in the structure of the Hofstadter spectrum is predicted to be undulatory 
particularly at the low value of t/ = 2t as shown in Fig. 3. For the sake fo clarity, undulatory 
feature of the band gap as a function of magnetic flux is explicitly shown in Fig. 5. It 
is noted that this undulatory behavior seen at a low value of C/ ~ 2t originates from the 
oscillatory variation of the staggered magnetization m as a function oip/q. At larger values 
of U such oscillatory behavior disappears as seen in Fig. 4. This is because the staggered 
magnetization (or antiferromagnetic order) becomes increasingly stable and thus insensitive 
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to the variation of magnetic field. The Hofstadter spectrum yields an additional sub-band 
only for odd q. This is caused by the Coulomb interaction which splits the band containing 
the fermi level into two bands. In short, one can readily see from the dispersion relation 
(17) that only the band gaps and the band widths in the Hofstadter spectrum but not the 
essential characters of the self-similarity and the Cantor set can be affected by interactions 
between electrons. 

3 Conclusion 

In the present study we derived a new generalized Harper's equation to deal with both the 
undoped (half-filling, 5 — Q) and doped {5 ^ 0) systems by allowing interactions between 
electrons, and examined the Hofstadter energy spectrum as a function of Coulomb repulsion 
U for the square lattice of correlated electrons at half-filling. We found that such Coulomb 
interactions between electrons affect only the band gaps and band widths in the structure of 
the Hofstadter butterfly. Judging from the dispersion relation (17), the inclusion of electron 
correlation effects beyond the mean field level will modify only the band gaps and band 
widths in the Hofstadter spectrum, but not the self-similarity of the Hofstadter spectrum 
and the homeomorphism to the Cantor set. 
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Figure Captions 

Fig. 1 Staggered Magnetization with and without magnetic field; sohd fine for the infinite 
size lattice and dots, for the finite size lattices. 

Fig. 2 Hofstadter spectrum of non-interacting electrons at half-filling. 

Fig. 3 Hofstadter spectrum of interacting electrons at half-filling with U — 2t. 

Fig. 4 Hofstadter spectrum of interacting electrons at half-filling with U = 4t. 

Fig. 5 Band gap A as a function of -. 
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